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Abstract 

Cloaking using complementary media was suggested by Lai et al. in m- This was 
proved in [21] in the quasistatic regime. One of the difficulties in the study of this problem 
is the appearance of the localized resonance, i.e., the fields blow up in some regions and 
remain bounded in some others as the loss goes to 0. To this end, the author introduced 
the technique of removing localized singularity and used a standard three spheres inequality. 

The method used in m also works for the Helmholtz equation. However, it requires small 
size of the cloaked region for large frequency due to the use of the (standard) three spheres 
inequality. In this paper, we give a proof of cloaking using complementary media in the 
finite frequency regime without imposing any condition on the cloaked region; hence the 
cloak works for all frequency. To successfully apply the approach in [21] , we establish a new 
three spheres inequality. A modification of the cloaking setting to obtain illusion optics is 
also discussed. 

MSC: 35B34, 35B35, 35B40, 35J05, 78A25, 78M35. 

Key words: cloaking, illusion optics, superlensing, three spheres inequality, localized reso¬ 
nance, negative index materials, complementary media. 

1 Introduction 

Negative index materials (NIMs) were investigated theoretically by Veselago in [36|. The 
existence of such materials was confirmed by Shelby, Smith, and Schultz in [35] . The study 
of NIMs has attracted a lot attention in the scientific community thanks to their interesting 
properties and applications. One of the appealing one is cloaking using complementary media. 

Cloaking using NIMs or more precisely cloaking using complementary media was suggested 
by Lai et al in m- Their work was inspired from the notion of complementary media suggested 
by Pendry and Ramakrishna in [32] . Cloaking using complementary media was established in [21] 
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in the quasistatic regime using slightly different schemes from m- Two difficulties in the study 
of cloaking using complementary media are as follows. Firstly, this problem is unstable since the 
equations describing the phenomenon have sign changing coefficients, hence the ellipticity is lost. 
Secondly, the localized resonance, i.e., the field blows up in some regions and remains bounded 
in some others, might appear. To handle these difficulties, in [21] the author introduced the 
technique of removing localized singularity and used a standard three spheres inequality. The 
approach in [21] also involved the reflecting technique introduced in [18]. The method in [21] 
also works for the Helmholtz equation; however since the largest radius in the (standard) three 
spheres inequality is small as frequency is large (see Section [2] for further discussion), the size of 
the cloaked region is required to be small for large frequency. 

In this paper, we present a proof of cloaking using complementary media in the finite fre¬ 
quency regime. Our goal is not to impose any condition on the size of the cloaked region 
(Theorem [T|) ; hence the cloak works for all frequency. To successfully apply the approach in 
|21j . we establish a new three spheres inequality for the second order elliptic equations which 
holds for arbitrary radius (Theorem [2] in Section [2]). This inequality is inspired from the unique 
continuation principle and its proof is in the spirit of Protter in [3l|. A modification of the 
cloaking setting to obtain illusion optics is discussed in in Section 0] (Theorem [3]) . This involves 
the idea of superlensing in m- Cloaking using complementary media for electromagnetic waves 
is investigated in [23] . 

In addition to cloaking using complementary media, other application of NIMs are super¬ 
lensing using complementary media as suggested in |29l |30l [33| (see also m) and confirmed 
in [niEsi, and cloaking via anomalous localized resonance m (see also laiioiEoi). Comple¬ 
mentary media were studied in a general setting in UHlEni and play an important role in these 
applications see naiiiiisniEiiEiiEsi. 

Let us describe the problem more precisely. Assume that the cloaked region is the annulus 
Byr 2 \ Br 2 for some r 2 > 0 and 1 < 7 < 2 in which the medium is characterized by a matrix 
a and a function a. The assumption on the cloaked region by all means imposes no restriction 
since any bounded set is a subset of such a region provided that the radius and the origin are 
appropriately chosen. The idea suggested by Lai et al. in m in two dimensions is to construct 
a complementary media in \ B^^ for some 0 < ri < r 2 . 

In this paper, instead of taking the schemes of Lai et al., we use a scheme from m which is 
inspired but different from the ones from m- Following [21], the cloak contains two parts. The 
first one, in Br 2 \ makes use of complementary media to cancel the effect of the cloaked 
region and the second one, in B ^^, is to fill the space which “disappears” from the cancellation by 
the homogeneous media. Concerning the first part, instead of B.yr 2 \ Br 2 , we consider Br^ \ Br 2 
with r 3 = 2 r 2 (the constant 2 considered here is just a matter of simple representation) as the 
cloaked region in which the medium is given by 

{ a, a mByr 2 \Br 2 , 

I, 1 in Br^ \ B^r 2 ■ 

The complementary media in \ B^ is given by 
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( 1 . 1 ) 


where F : 3^2 \ —>■ \ By .2 is the Kelvin transform with respect to 83^2, i-e., 

r 2 

F{x) = -r^x. 

fI 

Here 


T*a{y) = 


DT{x)a{x)DT^{x) 

J{x) 


and T^cr{y) = 


J{xy 


where x = T ^{y) and J{x) = | det-DT(x)| for a diffeomorphism T. It follows that 


n = 

Concerning the second part, the medium in B^^ is given by 

'fz!rt) h,(r=/ri)l 


The reason for this choice will be explained later. 

With the loss, the medium is characterized by sqS, where 


= < 


a, o 


F^ ^a, F^ ^(T 

rl/rlY h, (ri/r|)‘ 
FI 


in Bj.^ \ 3^2 , 
in Bj ~2 \ B^i , 


in B. 


ri) 


otherwise, 


and 


•Si = 


— l + i5 in Br 2 \ B ^, 
1 otherwise. 


( 1 . 2 ) 


(1.3) 


(1.4) 


(1.5) 


Physically, the imaginary part of is the loss of the medium (more precisely the loss of the 
medium in 3^2 \ 3^^). Here and in what follows, we assume that. 


< a{x)Y, ■ i VC £ for a.e. x G B^r 2 \ Fr 2 ^ 


( 1 . 6 ) 


for some A > 1. In what follows, we assume in addition that 


a is Lipschitz inH^-g \ (1-7) 

We can verify that medium sq^ is of reflecting complementary property, a concept introduced 
in |18l Definition 1], by considering diffeomorphism G : W^\Br^ B^^ \{0} which is the Kelvin 
transform with respect to dB^^, i.e.. 


G{x) = r|x/|xp. 


( 1 . 8 ) 


It is important to note that 


G^F^,A = / in Br^ 


(1.9) 
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since G o F{x) = This is the reason for choosing A in (II.3p . 

Let O be a smooth open subset of (d = 2, 3) such that CC ri. Given / G L^(0), let 
ri 5 , u G Hq{Q) be respectively the unique solution to 

div(s5^Vri5) + sok'^TiUs = f m Q, (1-10) 

and 

Au + k'^u = / in fl. (l-H) 

As in m, we assume that 

equation (|l.lll) with / = 0 ,has only zero solution mHQ{Q). ( 1 - 12 ) 

Our result on cloaking using complementary media is: 

Theorem 1. Let d = 2, 3, / G L‘^{Lt) with supp f C Ll\ and let u and us in Hq[Q) be the 

unique solution to (|1.10l) and (|l.lll) resp. There exists 70 > 1, depending only on A and the 
Lipschitz constant of d such that i/ 1 < 7 < 70 then 

Us ^ u weakly in \ Br^) as 6 ^ 0. (1-13) 

For an observer outside Br^, the medium in Bj.^ looks like the homogeneous one by (I1.13P 
(and also (jl.lip l: one has cloaking. 

Remark 1. The case k = 0 was established in |21j . 

The proof of Theorem [1] is given in Section [3l It is based on the removing localized singularity 
technique introduced in |21] and uses a new three sphere inequality (Theorem [2]) discussed in 
the next section. The discussion on illusion optics is given in Section 01 

2 Three spheres inequalities 

Let V be an holomorphic function defined in Bji^, Hadamard in [8] proved the following 
famous three spheres inequality: 

< h\\L<^(dBii^)h\\L°^(dBr^) ( 2 - 1 ) 

for all 0 < iii < i 22 < -^3 where 



A three spheres inequality for general elliptic equations was proved by Landis m using Carleman 
type estimates. Landis proved m Theorem 2.1] that 0 if u is a solution to 

div(MVu) + b ■ Vv + cu = 0 in Br^, (2.2) 

^In fact, | 13l Theorem 2.1] deals with the non-divergent form; however since M is assumed C^, the two forms 
are equivalent. 
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where M is elliptic, symmetric, and of class 6, c € C^, and c < 0, then there is a constant 
C > 0 such that 

(2-3) 

for some a € (0,1) depending only on R 2 /Ri, R 2 /R 3 , the ellipticity constant of M, and the 
regularity constants of M, b, and c. The assumption c < 0 is crucial and this is discussed in 
the next paragraph. Another proof was obtained by Agmon [1] in which he used the logarithmic 
convexity. Garofalo and Lin in [6] established similar results where the L°°-norm is replaced by 
the L^-norm, and M is of class C^, b and c are in L°°\ 

hWL^idBa^) ^ Ch\\h(dBii^)\\'^\\L'2{dBr^) 
using the frequency function. 

A typical example of (12.2p when c > 0 is the Helmholtz equation: 

Av + k‘^v = 0 in (2.5) 


Given k> 0 , neither (12.41) nor (12.3|) holds for all Ri < R 2 < Rs- Indeed, first consider the 
case d = 2. It is clear that for n G Z \ {0}, the function Jn{kr)e^'^^ is a solution to (j2.5l) in 
\ {0}, where is the Bessel function of order n. By taking i2i, R 2 ., and i ?3 such that 
Jn{kRi) = 0 / Jn{kR2), one reaches the fact that neither (j2.4p nor (12.3p is valid. The same 
conclusion holds in the higher dimensional case by similar arguments. In the case c > 0, (j2.4p 
holds under the smallness of R^ (see e.g., [21 Theorem 4.1]); this condition is equivalent to 
the smallness of c for a fixed R 3 by a scaling argument. 

In this paper, we establish a new type of three spheres inequalities without imposing the 
smallness condition on R 3 . This inequality will play an important role in the proof of TheoremUl 
Define 

IbllH(aBr) = \\''^\\m/^{dBr) + \\MVv ■ i^Wn-i/^idBr)- (2-6) 

Here and in what follows, z/ denotes the outward normal vector on a sphere. 

Our result on three spheres inequalities is: 


Theorem 2. Let d > 2, ci,C 2 > 0, 0 < < i22 < -^3 < R*, cind let M be a Lipschitz 

uniformly elliptie symmetrie matrix-valued funetion defined in Bji*. Assume v € H^{Bii^ \Bji-^) 
satisfies 

|div(MVu)| < ci|Vu| + C 2 |u|, in Bji^\ Bji-^. (2-7) 


Then, for any Aq > 1 with R 2 G (XqRi, R 3 /Xq), there exist a constant C and q >1, depending 
on the elliptic and the Lipsehitz constant of M, C also depends on ci, C 2 , and Xq sueh 


that 


'^WnidBR^) < C|l^llH(asHj)ll'^llH(aBfl3) 


( 2 . 8 ) 


where 


a := 


fCl fCg 


(2.9) 
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In Theorem [2l one does not impose any smallness condition on Ri, i?2, -R3 and the exponent 
a is independent of c\ and C2. The proof of Theorem [2] is inspired from the approach of Protter 
in |34] . Nevertheless, different test functions are used. The ones in [3Tj are too concentrated at 
0 and not suitable for our purpose. The connection between three spheres inequalities and the 
unique continuation principle, and the application of three spheres inequalities for the stability 
of Cauchy problems can be found in [2]. 

The proof of Theorem [2] is presented in the next two subsections. 


2.1 Preliminaries 


This section contains several lemmas used in the proof of Theorem [2j These lemmas are 
in the spirit of [34]. Nevertheless, the test functions used here are different from there. Let 
0 < i?i < i?3 < +00. In this section, we assume that M is a Lipschitz symmetric matrix-valued 
function defined in \ and satisfies 

j|eP<M(x)^e<A|^|2 


A' 

for a.e. x G Br^ \ Br^, for some A > 1. Set 

L:= IlMlIioc+i?3||VM||Loo. 

The first lemma is: 

Lemma 1. Let d> 2 and z G H^{Bji^ \ Br-^). We have 


( 2 . 10 ) 


[ _ (x- MVz) div(MVz) >- j 

J Bur,\B u. J B 


for some positive constant C depending only on d. 
Proof. An integration by parts gives 


CL^\Vz\^- / _ CLVlVzp, 

Jd{BR^\BR^) 


_ {x ■ MVz) dw{MVz) = — / _ V{x-MVz)-MVz 

Bn, JBii^XBji-^ 




Using the symmetry of M, we haved 


+ / _ {x ■ MVz) MVz ■ u. (2.11) 

Jd(BR^\BR^) 


A 

dxi 


d 

{x ■ MVz) = 7^ ( MkjXj 


and 


/ _ 2xjMkj 

Jbr^\Br^ 


dxi 

92 


dz 

dxk 


— -MkjXj 


92 ^ 


dz dMkj dz 


dxidxk 


+ Affci 7^-1“ X 

dxk 


dxi dxk 


( 2 . 12 ) 


•2; , ^ dz 

Mil— = - 


L 


dxidxk dxi 

d{xjMkjMii) dz dz 


d f dz dz\ 
dxi) 


f nT d f dz 

/ XjMkjMil— — 

Jbr^\Br^ dxk \dxi 

L 


Br,\Br^ dXk dXidxi Jd{BR^\BR^) 


dz dz 

_ XjMkjMil ——Uk. (2.13) 


' dxi dxi 


^In what follows, the repeated summation is used. 
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We derive from ()2.12p and (I2.13P that 

- [ _ \/{x- MVz) ■ MVz >- [ _ - [ 


The conclusion now follows from (12.111) and (j2.14l) . 
The second lemma is 


d{BR^\BR^) 


CL^r\Vz\\ (2.14) 


□ 


Lemma 2. Let d > 2, /3 G M, and z G \ Br^). There exists p\^L > 1 such that if 

P > PA,L o,nd \ld\Kf^ > 2 then 

[ _ V|z|2)div(MVe“^^”'’)+ [ _ CLV/^V-^P-ij^p 

f _ jA-V/32r-2p-2|2|2 - f _ CL‘^\Vzf. 

JbrABr, 2 Jbr^\Br^ 


> 


'Br^\Bji^ 

for some positive constant C depending only on d. 

Proof. A computation yields 

div(MVe“^'’ ’’) = pl3e~^'" ^ — (p + 2)r“P“^]x • Mx ^ div(Mx). 

An integration by parts gives 

ef^^~\Mx ■ V\zf) div{MVe-^^~'’) = P + Q. 


' Br^\Br-^ 


(2.15) 


Here 

with 


P = Pi + P 2 + P 3 


Pi 

P 2 

P 3 


- f _ p‘^P‘^\z\‘^dW[r-‘^P-^(x-Mx)Mx], 

J Br^\Br^ 

[ _ p/3(p + 2)|z|^div[r“^“^(a; • Mx)Mx], 

JBr^\Br^ 

f _ 2pf3r~P~‘^div{Mx)zVz ■ Mx], 

J Bro \-Br, 


and 


Q = 




pI3\z\'^ (^[p/3r —(p + 2)r p ^]x ■ Mx'^Mx ■ u. 


We next estimate P and Q. A computation yields 

—div[r“^P“'^(x • Mx)Mx] = {2p + 4)(x • Mx)‘^r~‘^P~^ — r“^^“^div[(x • Mx)Mx]. 
This implies 


Pi > 


[ _ p^p‘^r-^P-^\z\^[{2p + 4)A-^-CL'^]. 

J Br^\Br^ 


(2.16) 
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Similarly, 


(2.17) 


P 2 >- _ {p + 2)p\f3\r P ‘^\z\‘^[{p + A)A “^+01“^]. 

A combination of (j2.16p and (I2.17P yields 

Pi + P2> [ _ A~^p^(3^r-^P~‘^\zf. (2.18) 

Here we used the fact that p > pa,l and > 2. On the other hand, using Cauchy’s 

inequality, we have 

IP3I < / _ p‘^P^r-^P~^L‘^\z\‘^+ [ _ CL^lVzp. 

It follows from (I2.18P that 

P> [ _ lp^P^A-^r-^P-^\z\‘^ - [ _ CL‘^\Vzf, (2.19) 

2 Jbji^xbr^ 

provided that p > 2A?‘L?‘. Since 


\Q\< [ _ 2A^p^/3^r VP- 

J d{Bji^\Bii-^) 

the conclusion follows. □ 

Using Lemmas [T] and [21 we can prove the following result. 

Lemma 3. Let d > 2, /3 S M, and v € \ Br^). There exists a positive constant pa^l > 1 

such that if p > pa,l and \ j3\R^^ > 2 then 


L 


^p+2^2l3r-P 

2p|/3| 


[div(MVu)]^+ [ 

Jb. 


CL^e^^'^ '’iVup 


+ [ _ CLVVr-2p-ie2^’''"|r;p+ [ _ 

Jd{BR^\Bni) Jd{BR^\BRj 


> 


[ _ iA-VVr-2p-2e2/3--"|up, 

Jbr^\Br^ 2 


for some positive constant C depending only on d. 

Proof. By considering the real part and the imaginary part of v separately, one might assume 
that V is real. Set 

2 = equivalently v = e~^'' ^z. 

Since div[M'V{gh)) = 2'Vh-MVg + hdiv{MVg)+gdiv{MVh) {M is symmetric), it follows that 


div(MVu) = 2/3pr-P-2g-/3r V • MVz + Viv(MVz) + zdiv(MVe"^’' '’). 
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Using the inequality (a + 6 + c)^ > 2a{b + c), we obtain 

i[div(MVu)]^ > 2\l3\pr-P-^e-^''~\x ■ MVz) (^e-^^”’’div(MVz) + zdiv(MVe-^''”'’)) . 
This implies 


^p+2g2/3r-P 


[ - ——rdiv(MVu)l^ > / 2(x • MVz)div(MVz) 

2p|/3| 


+ / _ "{Mx ■V\z\'^)diY{MVe-^^ "). 


Applying Lemmas [T] and [2l we have 

r ^p+2 20r-P j- 

/ ——rdiv(MVu)l > / 

Jbr^\br^ 2p|/3| 

-L 


Br^\Br-^ 


- CL^\Vzf^ 


d{BR^\BR^) 


CL‘^p^P‘^r-^P-^\z\‘^ + CLV| . (2.20) 


Sinrt^ r = ^ 1) 

|Vz|2 < 2e2/'’'-"(|Vu|2 +p 2 ; 32 ^- 2 p- 2 |^| 2 ^_ 

A combination of (j2.2np and (I2.2ip yields, since p > pa,l, 


( 2 . 21 ) 


L 


^p+2g2/3r-P 


Br^\Br-^ 2p|/3| 


The conclusion follows. 
We also have 


[div(A/Vu)]^ > / _ -CL‘^\VvA 

Jbr^\Br^ V2 J 


ld{BR^\BR^) 


(ULV/3V-2p-1|u| 2 + CL^rlVul^) . 


□ 


Lemma 4. Let d > 2, /3 G M, and v G \ Br^). There exists a positive constant pa^l > 1 

such that if p > pa,l and \ j3\KfP > 2 then 


L 


'’u div(MVu) + / _ 

Br^Bri JBr^\Br^ 

<[ _ / _ Ce2/^’'“"(r|Vu|2 + r->|2), 

Jbr^\Br^ Ja{BR^\BR^) 


for some positive constant C depending only on d, A, and L. 
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Proof. We have 




div{M\/v) = [ _ MSIv ■ '’v)- [ _ e^^’' '’vMVv ■ v. 

J Br^XBr, d(Bfi^\Bii^) 

( 2 . 22 ) 


Jbr^XBii-^ JBr^XBr^ 

On the other hand, 

MVv ■ = [ 

Je 


J Br^XBr^ 

and 

Since 


'BrABr. 


L 


d{BR^\BR^) 


'’MV?;-V?;-2/3pr“P“2g2/3r \mVvxJ (2.23) 
- v < f e2^'’”'’(r|V?;|2 + LV->|2). (2.24) 


ld(BR^\BR^) 


1 


2l3pr-P-'^vMVv ■ x < -A“^|V?;p + 


we derive from (j2.22p . (I2.23D . and (12.241) that 

[ _ e^^''~'’vdiv{M\/v)+ [ _ 

JBr^\Br^ JBr^\Br^ 2 

<[ _ C/3Vr-^P-\‘^P^-^\v\^+ [ _ Ce^^^-^{r\Vv\^ + r-^\v\^). 

Jbr^\Br^ Jd{BR^\BR^) 

The conclusion follows. 

Combining the inequalities of Lemmas [3] and [H we obtain 


□ 


Lemma 5. Let d > 2, /3 G M, and v G \ There exists a positive constant p\^L > 1 

such that if p > pa,l and \j3\ > 2Rf^ then 

[ _ + <C ^ \d\y{MVv)\^ 

Jbr^\Br^ ^ ^ JbrABr, 


Br^Bri 

+ C [ _ |/3|pe2^''”'’(r|V?;p+p2/?V-2p-i|?;|2). (2.25) 

Jd(BR^\BR^) 


for some positive constant C depending only on d, A, and L. 

Proof. Note that 

|?;div(MV?;)| < p|/3|>|div(MV?;)|VP+2. 

P\P\ 

The conclusion now follows from Lemmas [3] and HI The details are left to the reader. □ 
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2.2 Proof of Theorem [2] 

Let 

1 < A < Aq, 

(which will be defined later) and set 


D - Bxrs \ 

Let ui € H^{D \ BBr^) and tts € H^{D \ BBr^) be respectively the unique solution to 


/ 

div(MV/ii) = 0 

in D \ BBr^ 

[ui] 

= v] [MVui • r] = MVu • R 

on OBr^ , 

V 

ui = 0 

on clD, 


div(MVtt3) = 0 

in D \ Br^ , 

^ [us] 

= v] [M'Vu 3 ■ v] = M'Vv ■ V 

on dBR ^, 

K 

CO 

II 

O 

on dD. 


Here and in what follows, [•] denotes the jump across a sphere and v denotes the unit outward 
normal vector on a sphere. It follows that 


'^l\\m{D\dBR^) < 11^111/^3/2(95^^^^) < C'lbllH(aBHj) (2.26) 


and 

IWsW H^D\dBii^) < C'lbllH(aSH3), 11^311 i/3/2 (95^^^) < C\\v\\iii^SBii^)- ( 2 - 27 ) 

Here and in what follows in this proof, C denotes a positive constant depending only on the 
elliptic and the Lipschitz constant of M, ci, C2, Aq, R*, R*, and d. Set 


di — (A — 1 ) 77 / 

Let ipi, € C'c(M'^) be such that 

1 RRi+di/3\RRi^ 

0 in \ {BBi+di /2 \ Br-^/x) 

Define 

f V - ^PlUl - ip3U3 mUR^\ljR^, 

V = < _ (2.28) 

[ -ipiUi-ip 3 U 3 in D\{Br^\Br^). 

Applying Lemma [5l we obtain, for l/3| > 2 ( 7773 )^ 

C [ e2^’'“'’/3(/32|H|2 + |VH|2) < [ e2^’'"’’|div(MVH)|2+ [ {\VV\‘^+P^\V\‘^). (2.29) 

Jd Jd JdD 

The proof is now quite standard and divided into two cases. 



and (73 = (A — 1)7?3/A. 


and ip 3 = 


1 in 77^3 \ 77^3 _^ 3 / 3 , 

0 in \ (77 a_R 3 \ 7753_rf3/2)- 
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Case 1 : ||i'||H(aBflj) < lbllH(( 9 Sfl 3 )- We deduce from ()2.29p that for /3 > /3o, 


C + |VF|2) < (2.30) 

where 

/?3 = R^/X and Ri = RifX. 

This implies 

C\\v\\u{dBu^) < (2.31) 

Define a' € (0,1) and /3 > 0 as follows 

a' = ^ and /3(i?“^ - R~^) = (1 - a') In (||^;||H{aBfl 3 )/ll^^llH{aBflj) I- 

-^1 “ 213 


Note that 0 < a' < 1 since R 2 < Rz/l- We assume that ||u||H(aBfl 3 ) > for some 

large C such that /? > max{2i2^^, 2,/3o} since if ||n||H( 9 Sfl 3 ) < C'l|r’||H(aSflj)) the conclusion 
holds for any a £ (0,1) by taking /3 = max{2i?)^^, 2,/3o} in (I2.3ip . It follows from p2.3ip and 
the choice of a' and j3 that 


lbllH(aBjj2) < C'/^ll^llH(aBfl^)ll^llH(aBfl3)- 

Define 

_ 

Rl R^ 

It is clear that a < Hence, by choosing A close to 1, 


(2.32) 

(2.33) 


a < a' yR2 € {jqRi, Rs/'jo)- 


(2.34) 


A combination of p2.32l) and (12.341) implies 


(2.35) 


Case 2 : ||n||H(aB^j) > l|r’||H(aBfl 3 ). The proof is similar to the previous case by considering 
/3 < —2(7i?3)“^’. The details are left to the reader. The proof is complete. □ 


3 Cloaking using complementary media. Proof of Theorem [T] 

This section containing three subsections is devoted to the proof of Theorem [H In the 
first subsection, we present two useful lemmas. The proof of Theorem [1] is given in the second 
subsection. 

^Here we assume that ||u||H(8Siij) A 0 since otherwise v = 0. This fact is a consequence of the unique 
continuation principle and can be obtained from (12.3111 by letting /3 —>■ oo. 
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3.1 Preliminaries 

In this section, we present two lemmas which will be used in the proof of Theorems [U and El 
The first lemma is on a change of variables and follows from |18l Lemma 1]. 

Lemma 6. Let d > 2, k > 0, 0 < Ri < R 2 < Rs with R^ = R 2 /R 1 , a G \ 

be a matrix valued function, a G L°°{Bii^ \ Br^) he a complex function, and K : Br^ \ Br^ —)• 

Br^ \ Br^ be the Kelvin transform with respect to OBr^, i.e., 

K{x) = R 2 x/\xf. 

For V G H^{Br^ \ Br^), define w = v o K~^. Then 

div(oVu) + k'^av = 0 in Br^ \ Br^ 

if and only if 

d\v{K^aVw) + k^K^fUw = 0 m Br^ \ Br^. 

Moreover, 

w = V and K^.a'Vw ■ v = —aSJv ■ v on OBr^. 

The second lemma is a stability estimate for solutions of (|1.10l) . 

Lemma 7. Let 0 < 5 < 1, f G L‘^{Ll), and let A G and S G L°°{Q,C) be such 

that A is Lipschitz and uniformly elliptic and ^(S) > 0. Assume that us G Hq(LI) is the unique 
solution to (ll.lOp . Then 

ll'^'5ll^i(n) — Cl (^S ^||/||l 2 (q)+ ||/|li2(Q)^ , (3-1) 

for some positive constant C independent of 6 and f. 

Lemma [7] is a variant of [181 Lemma 1]. The case A: = 0 and its variant in the case k > 0 
were considered in [21] and m respectively. The proof is similar to the one of [181 Lemma 1]. 
For the convenience of the reader, we present the proof. 

Proof. The existence and uniqueness of us are given in |18j . We only establish (EH) by 
contradiction. Assume that dMD is not true. Then there exist —>■ 0 and (/„) C L?‘{Q) such 
that ^ 

ll'“n||_H-i(n) = 1 and T-\\fn\\L'^[n)\Wn\\L'^{suppfn) + \\fn\\\ 2 (a) ^ 0, (3.2) 

as n —>■ oo, where Un G Hq{Li) is the unique solution to 

div(s5^AV'u„) + k‘^soT,Un = fn in P. (3.3) 

Without loss of generality, one may assume that Un ^ u weakly in H^{Q) and strongly in Lp‘{Ll)\ 
moreover, u G Hq{LI) and u satishes 

div(soAVu) + k^soTjU = 0 in P. (3-4) 
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Multiplying equation ()3.3p by Un (the conjugate of Un) and integrating on n, we have 

/ Ss^AVUn-VUndx - / SQ'Z.\Un\^ dx = - / fnUndx. 

Jvt JQ. JQ. 

Considering the imaginary part and using the fact that 


1 


/ 

Jn 


fnUr 


1 

< - 


^„.,JlL2(n)lkn||L2(snpp/„) ^ 0 as U CX) by (I32|), 

liVu,; 


0 as n —>■ oo. 


(3.5) 


we obtain, by (II.6p . 

Since div(AVun) + k'^TiUn = fn in Br 2 \ Br^ and /n ^ 0 in L^(n), it follows from that Un ^ 0 
in the distributional sense. This in turn implies 


hnh^Br^XBri) ^ 0 as n ^ oo. 
A combination of (j3.5p and (j3.6p yields 

hnWmiBr^XBr^) ^ 0 as n ^ oo. 

Hence 


(3.6) 

(3.7) 


u = 0 in Br 2 \ B, 


ri ) 


and 


(3.8) 

Since u = 0 in \ Br^ and u satisfies (|3.4p , it follows from the unique continuation principle 
that tt = 0 in H. Hence, since —>■ u in L^(0), 


Un ^ 0 in L^(17) as n —>■ oo. 

Multiplying (|3.3p by Un and integrating on 0 \ Br 2 , we have 

/ AVUn'VUndx— / k'^SoT.\Un\'^ dx = — / fnUndx + 

J Q,\Br2 Q\Br2 ^ 

Using (13.8p and (j3.9p . we obtain 

||Vun||L 2 (Q\s^j ^ 0 as n ^ OO. 

Similarly, 

||Vun||i, 2 (g^^) 0 as n oo. 

A combination of (j3.7p . (|3.9p . (I3.10p . and (13.lip , we obtain 

ll^nll^qn) —>■ 0 as n ^ oo. 
which contradicts ()3.2p . The proof is complete. 


(3.9) 


AVu 


n ■ ^2 Un- 


'dBr 


(3.10) 

(3.11) 


□ 


14 








3.2 Proof of Theorem [T] 

We use the approach in [2T] with some modifications from m so that the same proof also 
give the result on illusion optics (Theorem [3] in Section [3|). However, instead of applying the 
standard three sphere inequality as in m , we use Theorem [21 
We have, by Lemma [71 


||«5||^i(o) < C \^6 ^\\f\\L‘ 2 {n)\\us\\L 2 (^n\Br^) + " 
As in [21], let ui^s be the reflection of us through dBr 2 by F, i.e. 

^^1,5 = Us o F~^ in \ Br 2 
and let M 2,5 be the reflection of ui^s through dBr 2 by G, i.e., 

U2,5 = 'Ui,<5 o G~^ in Br^. 

By Lemma [6l 


1 


div(AVui^ 5 ) + -- —k^TiUi^s = 0 in Br^ \ B, 


1 — i5 

u 2 


Au2,5 + k U 2 ^s = 0 in B^^. 


r 2 ) 


■'rs \ 1 


Ul,S = Us 


on dBr 2 and (1 — i6)A'Vui^s ' ^ = AVus ■ v 


1 

(3.12) 


(3.13) 


(3.14) 


(3.15) 

we have 

(3.16) 

on dBr 2 - 

(3.17) 


Let Vi^s S \ Brj) be the unique solution to 

iS 

div(AVVi_5) + k'^'FVi^s = - r^k'^'Bui^s in \ Hj-j, 

• u — ikVi^s = 0 on dBr 2 , 

Pi,5 = 0 

By Fredholm’s theory, 

l|Pl,<5||i^i(W3\W2) - 11'“<5II Hi (!^)- 

Define Ui^s ia B^^ \ follows 

Ui^s = us — Ui^s — Pl,(5- 


on dBr^. 


(3.18) 

(3.19) 

(3.20) 


Then Ui^s S H^{Br^ \ H^j) and Ui^s satisfies 

div(AV17i^5) + k^TiUi^s = 0 in Br^ \ Hj-j 
l|Pi,<5||_H-i/2(aB^^) + \\AVUi^s ■ ^ll_ffi/2(a_B^2) - ^'^ll“<5llHi(n)! 


and 

\\Bl,s\\H^/^(dBr^) + ll^^f^l,<5 • ^ll_ffi/2(a_B^g) < <71111511^^1(^2). 
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Applying Theorem [2l we have 

\Wl,5\\H^/2[dB-^r2) ^ ■ '^11^1/2(95^^^ < C5°‘\\us\\Hi{n) (3-21) 

where a is given in ()2.9p with Ri = r 2 , R 2 = 75 " 2 , R 3 = r^. By first taking q large enough and 
then choosing 70 close to 1 enough, from (j2.9p . we can assume that 

a > 1/2. (3.22) 

Here is the place where the condition 7 < 70 is required. A combination of (13.191) and (13.211) 
yields 


||^^<5 - Ul,s\\Hi-/2(dBjr2) “*■ ll^^(^<5 “ "“ 1 , 5 ) ■ ^\\H-^/^{dBjr2) - 11'“<511 Hi (!^)- (3.23) 

In what follows, we assume that k = 1 for notational ease. Define U 2 ,s in \ B^r 2 ns follows 


U 2 ,s = ui,s - U 2 ,s + yi, 5 - 


Then 


A[/ 2,5 + = 0 in Br^ \ B-yr 2 , (3.24) 

and 

tS 

1/25 = 0 and drU 2 s = —- r^drUi 5 + drVi 5 on dB^. (3.25) 

’ ’ 1 — zo 

Case 1 : d = 2 . As in m, define 

^ ^ TTZ 

Jnir) = 2”n!Jn(r) and Ynir) = _ ^y ^nir), 


where Jn and are the Bessel and Neumann functions of order n. It follows from m (3-57) 
and (3.58)] that 


Jn{t) = t”[l + 0 ( 1 )] 

(3.26) 

and 


Ynit) = t "'[1 + 0 ( 1 )], 

(3.27) 

as n —t + 00 . From (|3.24p one can represent Uo 5 as follows 



7 / 2,5 = aoJo{\x\) + 6ok"o(|a;|) + EE [an,±Jni\x\) + bn,±Yn{\x\)]e^^'^^ in Br^ \ B-yr^, (3.28) 

n=l it 

for ao,bo,an,±,bn^± G C (n > 1). Assume that 

drU 2 ,s = Co + ^ ^ on 55^3. 

n>l ib 
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Then, by (|3.18p . (|3.19p . and (j3.25p . 


+ ^|Cn,±P ~ ||5,C/2,5||^_l/2(gB,3) 

n>l di 

Using p3.25p again, we have 

f ®n,±dn(?’3) T ^n,±^^(^3) — 0 

< ^ ^ for n > 0. 

I an,±J'n{rz) + bn,±Y^{r3) = Cn,± 

Here we denote ao,± = oo/2; &o,± = &o/2; and co,± = co/2. It follows that 

I ®n,± — 


(3.29) 


^n,± — Cn^±BC,fi^ 


for n > 0. 


(3.30) 


where 


AC„ = — 


U, 


TV— T'Y 


(rs) and BCn = - 


Jn 


Y T' —Y' T 


{fz)- 


Using (I3.26P and (|3.27l) . we derive that 

ACn, = + °(^)) 

We now make use the removing of localized singularity technique introduced in |19l [2T] . Set 


^ ^6„,±y„(|x|)e=^*^® in Br^ \ B^r 2 , 

n=l it 

We claim that, for yr 2 < r < r^, 

\\U2,S — "“<511^1/2(95^) + \\drU2,5 “ '9rW<5|l5-i/2(95^) < CS\\us\\Hi-{Q)- 
Indeed, for 'yr 2 < r < r^, 


(3.31) 


\\U 2,5 "“< 5115 - 1 / 2 ( 95 ^) —II ^ ^ “n,±dn 


\x\]e 


inO II2 


lHl/2(9Br) 


+ l)l“ii,±l |dn(|* 


+ l)\cn,±ACn\^\Jn{\x\)? < C”X]Z] ^\cn,±?{r/■ 

n.>0 di 7T.>0 di 


It follows from (I3.29P that 


for yr 2 < r < r^. Similarly, 


||t^2,,5 - "“5|l5i/2(9BU ^ C'<^II"“<5 ||hi(0), 


\\drU2,S — 5r'“<5||5-i/2(95^) < C'<5||//(5||5i(f1); 
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for 7 r 2 < r < r^. As a consequence of (|3.19p and (I3.3ip . we obtain for 'yr2 < r < r^, 
hi,5 - U2,5 - Us\\Hi/2{dBr) + \\drUi,5 - du2,S “ ^rU5\\H-l/2^^QBr) ^ C6\\us\\m{n)- 
Define 


Us 


us in D \ Br ^, 

< Us — Us if X G Br^ \ B^r2 ) 
U2,S if X G Byr2 ■ 


We have 

div{AVUs) + k'^^Us = / in D \ {dBr^ U dB^^h- 
On the other hand, from (I3.23P and (I3.32p . we obtain 


(3.32) 


+ Wh^Us ■ < CdhusWmin) 

and 

Using p3.12p . we derive that 

Since a > 1/2, it follows that Us is bounded in \ {dBr^ U dB^r 2 ))- Without loss of 

generality, one may assume that Us ^ U weakly in {dBr^ UdB-yr 2 )) as 5 —>■ 0; moreover, 

U G H^{n) and 

AU + k'^U = / in D and U = 0 on (90. 

Hence U = u. Since the limit is unique, we have the convergence for the family (Us) as (5 —>■ 0. 
Case 2: d = 3. Define 


Jn{t) = 1 • 3 • • • (2n + l)jnit) and ijn = ^ 

1 • 3 • • • (2n — 1) 

where jn and yn are the spherical Bessel and Neumann functions of order n. Then, as n large 
enough, (see, e.g., [5l (2.37) and (2.38)]) 

jnikr) = r'^[l + 0{l/n)) and yn{kr) = + 0{l/n)). (3.33) 

Thus one can represent U 2 ,s of the form 

oo n 

^2,5 = (3.34) 

n=l —n 

for G C and x = x/|x|. The proof now follows similarly as in the case d = 2. The details 

are left to the reader. □ 
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Remark 2. In the proof, we use essentially the fact {A, S) = (/, 1) in \B-yr 2 to use separation 
of variables in this region. In fact, this condition is not necessary by using the technique of 
separation of variables for a general structure in [20]. 

Remark 3. The construction of the cloak given by ()1.4p is not restricted to the Kelvin trans¬ 
forms F (and G). In fact, one can extend this construction to a general class of reflections 
considered in m- 

Remark 4. The condition {F^,A, iTS) = [A, S) in \ Br 2 is necessary to ensure that cloaking 
can be achieved and the localized resonance might take place see |2l| (see also |1| for related 
results). 

Remark 5. Cloaking can also be achieved via schemes generated by changes of variables (TJ 
laEi]. Resonance might also appear in this context but for specific frequencies see UM- It 
is shown in [T6| that in the resonance case cloaking might not be achieved and the field inside 
the cloaked region can depend on the field outside. Cloaking can also be achieved in the time 
regime via change of variables [261 [27] . 


4 Illusion optics using complementary media 


We next discuss briefly how to obtain illusion optics in the spirit of Lai et al. in [12]. The 
scheme used here is a combination of the ones used for cloaking and superlensing in [211ES] and 
is slightly different from [12]. More precisely, set 


2 / 2 

= r3/r2. 

Let Uc G 1*6 elliptic and Uc G L°°(R^ 2 ^^ 2 ,C) with 9(crc) > 0. Define 

f A, S in 0 \ 


^ 1 , Si — 


ill jyn 1 


(4.1) 


and 


4i,Si = < 

\ {n/r2f 

Recall that (4, S) is defined in (11.41) . 
solution in Rq(0): 


/, 1 in D \ Bro, 

(4.2) 

ac{x/m),{r^/r2) ^Ocixlm) in Br^- 
We assume that the following equation has only zero 


div(^iVu) -|- fe^Siu = 0 in 0. (4.3) 

We obtain the following result on illusion optics: 

Theorem 3. Let d = 2, 3, / G L^(D) with supp / C D \ By^ and let u and us in Rg(D) he 
respectively the unique solution of 

div(s54iVu5) + k‘^soT,ius = / in Ll, 


and 

div(^i Vu) + fc^Siu = / inLl. 
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There exists 70 > 1, depending only on A and the Lipschitz constant of a such that ifl < 7 < 70 
then 

Us ^ u weakly in \ Bj.^) as (5 —)■ 0. (4-4) 

For an observer outside the medium in B^-^ looks like (Ai,I]i): one has illusion optics. 

Proof. The proof is similar to the one of Theorem [TJ Note that in the proof of Theorem [H we 
do not use the information of the medium inside The details are left to the reader. □ 
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